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ABSTRACT 

Aims. We investigate metallicity and a-element abundance gradients along a Galactic longitude strip, at latitude b-4°, with the 

aim of providing observational constraints for the structure and origin of the Milky Way bulge. 

Methods. High resolution (R~22,500) spectra for 400 K giants, in four fields within -4.8° < b < -3.4° and -10° < 1 < +10°, were 
obtained within the GIRAFFE Inner Bulge Survey (GIBS) project. To this sample we added another ~ 400 stars in Baade’s Window at 
(1, b) = (1°, -4°), observed with the identical instrumental configuration: FLAMES GIRAFFE in Medusa mode with HR13 setup. All 
target stars lie within the red clump of the bulge color magnitude diagram, thus minimizing contamination from the disc or halo stars. 
The spectroscopic stellar surface parameters were derived with an automatic method based on the GALA code, while the [Ca/Fe] 
and [Mg/Fe] abundances as a function of [Fe/H] were derived through a comparison with the synthetic spectra using MOOG. We 
constructed the metallicity distributions for the entire sample, as well as for each field individually, in order to investigate the presence 
of gradients or field-to-field variations in the shape of the distributions. 

Results. The metallicity distributions in the five fields are consistent with being drawn from a single parent population, indicating the 
absence of a gradient along the major axis of the Galactic bar. The global metallicity distribution is well fitted by two Gaussians. The 
metal poor component is rather broad, with a mean at < [Fe/H] >= -0.31 dex and <x = 0.31 dex. The metal-rich one is narrower, with 
mean < [Fe/H] >= +0.26 and cr = 0.2 dex. The [Mg/Fe] ratio follows a tight trend with [Fe/H], with enhancement with respect to 
solar in the metal-poor regime, similar to the one observed for giant stars in the local thick disc. [Ca/Fe] abundances follow a similar 
trend, but with a considerably larger scatter than [Mg/Fe]. A decrease in [Mg/Fe] is observed at [Fe/H] = -0.44 dex. This knee is in 
agreement with our previous bulge study of K-giants along the minor axis, but is 0.1 dex lower in metallicity than the one reported 
for the Bulge microlensed dwarf and sub-giant stars. We found no variation in tr-element abundance distributions between different 
fields. 

Key words. Galaxy: bulge - Galaxy: abundances - Galaxy: structure - Galaxy: evolution - Galaxy: formation - 


1. Introduction 


In the last few years we have witnessed important progress to¬ 
wards the understanding of the Milky Way bulge. The availabil¬ 
ity of large photometric and spectroscopic datasets from ongo¬ 
ing surveys are finally allowing us to obtain a wider view of the 
Bulge properties, even expanding our knowledge to the inner¬ 
most regions previously not accessible. 

The bar in the inner regions of the Galaxy was first suggested 
by de Vaucouleurs (1964) and its main structural properties have 
been investigated in detail since then ( Blitz & Spergel||1991| 
IStanek et al.|| 19941 |Dwek et al.|| 19951 IBabusiaux & Gilmorel 


* Based on observations taken with ESO telescopes at the La Silla 
Paranal Observatory under programme IDs 187.B-909(A) and 71.B- 
0196. 


|2005[ |Rattenbury et al.||2007| Gonzalez et al.||2~011| and ref¬ 
erences therein). The use of red clump (RC) giant stars as a dis¬ 
tance indicator has been a fundamental tool for Bulge morphol¬ 
ogy studies (iStanek et al.j 1994) . Based on this technique the ax¬ 
ial ratios of the bar have been currently constrained to be about 
1:0.4:03 with a bar size of about 3.1-3.5 kpc major-axis length. 
The position angle of the bar has been historically measured to 
range between ~20-40 deg. with respect to the Sun-centre line 
of sight with its near end towards positive Galactic longitudes 


dBlitz & Spergel 1991; Stanek et al. 1994;|Dwek et al.|[l995; 


Binney et al.||i997f Bissantz & Gerhard 2002; Benjamin et al. 


2005 Babusiaux & Gilmore|2005j Rattenbury et al.|2007 Cao 


et al. 201 3[ >. Recently, detailed comparison of RC stellar counts 

to N-body models (e.g. Wegg & Gerhard 2013) have allowed 
to constrain the bar position angle to 27 - 33°. The RC tech- 
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nique to measure distances towards the Bulge led to the discov¬ 
ery of a split RC at 1 = 0° and |b| > 5° by Nataf et al. ( |2010| > 


and Me William & Zoccali (2010 1 . This discovery was quickly 
followed by the construction of wider, more detailed 3D maps 
that revealed the X-shaped structure of the Bulge (Saito et al. 
2011 1 Wegg & Gerhard|20l3| ) as first suggested by |McWilliam 


& Zoccali|(|2010]>. This X-shape structure is often seen in exter¬ 


nal galaxies and is well reproduced by dynamical models of disc 
galaxies (e.g. |Li & Shen|20l2} Gardner et al.|2014 1 . They corre¬ 
spond to an extreme case of the boxy/peanut (B/P) structures 
observed in several external disc galaxies formed as the con¬ 
sequence of buckling instabilities of galaxy bars. The buckling 
process of the bar results in the heating up of stellar orbits in the 
vertical direction (e.g.,|Combes et al.|1990| Athanassoula 2005 1 


Mart inez-Valpuesta et al.||2006[ |Debattista et al.||2005| >. Due to 

the absence of an external influence in this process, a reference 
of these structures as so-called pseudo-bulges (see Kormendy & 
|Kennicutt|2004| for a detail definition of such structures) is of¬ 
ten found in the literature. However, we refrain from using such 
a definition to differentiate them from the young stellar disc-like 
structures formed in the inner part of bars due to the continuous 
in-fall of gas. We will then refer to this structure simply as the 
B/P bulge of the Milky Way. 


The Bulge Radial Velocity Assay survey (BRAVA; Howard 
et al.||2009l |Kunder et al.||2012) , measured radial velocities for 
- 10,000 Bulge M giants at latitudes b = -4°, -6°, -8° and lon¬ 
gitudes -10° <1 < 10°. The BRAVA survey provided strong ev¬ 
idence for the cylindrical rotation of the Bulge, concluding that 
a B/P bulge would be enough to reproduce the rotation curve of 
the Bulge without the need of a classical component, formed via 
mergers in the early evolution of the Galaxy (Howard et al.|2009 1 


Shen et al. 


2010] ). Moreover, adaptive optics resolved Ha kinematic studies 
have revealed that most z~2 main-sequence galaxies are large 
rotation-dominated discs, with a minor fraction of major merg¬ 
ers, are characterised by very high velocity dispersion (-50-100 
km/s) and by the presence of several kpc scale, actively star- 
forming clumps (e.g. [Forster SchreibeFet al. 2009; Genzel et al. 


2011, Mancini et al. |2011 [ Newman et a l.|2013) >.~In parallel with 

these findings, a new paradigm has emerged in which bulges 
form by the migration and coalescence to the center of massive 
star-forming clumps (jlmmeli et al.|2004[|Bournaud et al.[2 007 ; 
|Elmegreen et al.|2008[|Boumaud e t al. 2009]), a situation which 

may even develop into a global violent disc instability, leading 


to a dissipational formation of compact bulges (Dekel & Burkert 
|2014] >. Considering all the process that can be involved in the for¬ 
mation of bulges, it becomes clear that kinematics alone would 
not be sufficient to reconstruct the formation mechanism of the 
Galactic bulge. Indeed, the chemical abundance measurements 


of individual stars have been shown to provide a different per¬ 
spective, especially when combined with the kinematics. 

Several studies have investigated the metallicity distribution 
functions near to Baade’s Window at (l,b) = (0°,-4°) (e.g. 


Rich 

50 

oo 

00 

McWilliam & Rich 

1994a 

Minniti 1996; Sadler 

et al. 

1996 

Ramirez et al. 2000; 

Zocca 

i et al. 2003 

Fulbright 

et al. 

2006 

deriving a wide distri 

rution with metallicities rang- 


ing from -1.6 to 0.5 dex and a peak at solar metallicities. Zoccali 
et al.| ( 2008] ), and |Johnson et ah ( |201 1 2013 1 used high res 


olution spectra to derive metallicity distributions in different 
Bulge regions along the minor axis and have firmly established 
the presence of a vertical metallicity gradient with field stars 
closer to the Galactic center being on average more metal-rich. 
Furthermore, the general picture of these metallicity gradients in 
the Bulge was presented in Gonzalez et al.| ( [2013[ > who showed 
the global mean photometric metallicity map of the Galactic 
bulge based on the Vista Variables in the Via Lactea (V V V) ESO 
public survey data ( Minniti et al.|2010] ). 

By combining [Fe/H] and [Mg/Fe] abundances and kinemat¬ 
ics, Hill et al. ( 201 1] ) and Babusiaux et al. (2010) suggested the 
presence of two distinct components in the Bulge: a metal poor 

one (< [Fe/H] >-0.3) with kinematics typical of a classical 

spheroid, and a metal rich one, (< [Fe/H] >~ +0.3) concen¬ 
trated towards the Galactic plane with a significant vertex devi¬ 
ation, suggestive of a bar-like component (see also Soto et al. 


_ 2010| ). The cylindrical rotation was also confirmed re¬ 
cently by the GIBS survey in |Zoccali et al. ( 2014| hereafter Paper 
I) where it was extended to latitudes (b-2°). 

While the B/P bulge has been shown to rotate cylindrically, 
the bulk of its stellar population is over -10 Gyr old (e.g.|Zoccali| 
[eFaTl[2003l jClarkson et al.l[20TT] [Valenti et al.|[20T3] >. A -10 
Gyr look-back time brings us to redshift z~2, where galaxies are 
radically different from local ones and live at the epoch when 
the overall star formation rate (SFR) in the Universe peaked. 
Massive galaxies with stellar masses comparable to those seen 
today (~ 10 10 - 10 n Msun) are forming stars at rates of -20- 
200 Msun/yr, some 20 times higher than in local galaxies of the 
same mass. These higher SFRs are a consequence of higher gas 
mass fractions in these systems, typically -30-50% as revealed 
by CO observations (Tacconi et al. |2010| 2013] Daddi et al. 


2007). The origin of the observed metallicity gradients would 
then be the natural consequence of the different contribution of 
each of these components as a function of Galactic latitude. A 
similar conclusion was reached by Rojas-Arriagada et al. ( 2014) ) 
based on the Gaia-ESO survey data. 

A step further in the interpretation of the Bulge metallicity 
distribution was done using the Abundances and Radial velocity 
Galactic Origins Survey (ARGOS; Freeman et al.|[20l3 i. Ness 


|etal.| ( |20T3 a b) measured radial velocities, [Fe/H], and [a/Fe] ra 
tios for -28,000 Bulge stars at different regions well-distributed 
across the Bulge. They find the same cylindrical rotation found 
by the BRAVA and GIBS survey, but suggested the presence 
of 5 components in the metallicity distribution of the Bulge. 
They suggested that the actual Bulge component would be the 
metal rich one ([Fe/H] - +0.15), and perhaps the component at 
[Fe/H] - -0.25 dominating at high latitudes. Additional com¬ 
ponents would belong to the inner disc and halo. 

On the other hand, a-element abundances in Bulge stars pro¬ 
vide us with with an additional constraint on the formation his¬ 
tory of the Bulge stellar populations: its formation time-scale. 
As suggested by Tinsley ( |1979| the ratio of [a/Fe] compared to 
[Fe/H] is expected to be a function of the time delay between 
the production of both a- and iron-peak elements by SNe II (e.g. 
Woosley & WeaveT p 995) and the yield of mostly iron-peak ele¬ 
ments with little a-element production by SNe la (e.g.|Nomoto 


et al. 1984). Clearly, a sufficient amount of time needs to be ac¬ 


counted for such that enough SNe la events occur for the [a/Fe] 
ratio to decline from the SNe II value. The a-element abun¬ 
dances of Bulge stars with [Fe/H]<-0.3 have been historically 
found to be enhanced over iron by [a/Fe]—1-0.3 dex (]McWilliam| 
& Rich||1994b[ [Rich & Origlia| [20051 |Cunha & Smith|)2006[ 


Fulbright et al. 2007 |Fecureur et al. 2007| |Rich et al. 2007 

thus suggesting a fast formation scenario. Metal-rich stars on 
the other hand show a decrease in [a/Fe] reaching [a/Fe]= 0 for 
metallicities above Solar values. However, the SNe la delay time 
might vary depending on the different production channels that 
could be present, and therefore to obtain a direct translation of 
these trends to absolute time scales is not straightforward. As 
a consequence, the comparison of [a/Fe] values in Bulge stars 


2 










































































































































































O. A. Gonzalez et al.: The Giraffe Inner Bulge Survey (GIBS) II 


Table 1 . Observed fields and final number of stars analysed in this work. 


Field name 

RA 

(hr) 

DEC 

(deg) 

1 

(deg) 

b 

(deg) 

Setup 

R=/l/A/l 

A coverage 

A 

N s tars 

GIRAFFE 

exptime/star 

(s) 

Ay 

mag 

HRp8m3 

18:16:40.8 

-23:45:32.20 

7.9460 

-3.4770 

HR13 

22500 

6120-6405 

103 

27000 

1.71 

HRp4m3 

18:07:15.4 

-27:31:21.70 

3.6174 

-3.4111 

” 

” 

” 

88 

27000 

1.61 

HRm5m3 

17:47:49.2 

-35:03:24.10 

355.0036 

-3.5701 

” 

” 

” 

103 

27000 

1.59 

HRm7m4 

17:48:11.0 

-37:09:25.30 

353.2336 

-4.7106 

” 

” 

” 

106 

27000 

1.33 


against those of other galactic components has been proven to be 
a useful tool to investigate the Bulge formation time-scale with 
respect to those other components such as disc and halo. Based 
on this relative measurement approach, Fulbright et al. ( |2007j i, 
|Zoccali et al.| ( |2006l >, and |Lecureur et al.| ( |2007| ) concluded that 
the [ff/Fe] ratio was enhanced by nearly +0.1 dex with respect to 
the trends of both the local thin and the thick disc as traced by 
nearby dwarf stars. Therefore, these results implied a shorter for¬ 
mation time scale for the Bulge than from both discs. However, 
|Melendez et al. ( 2008 1 , |Alves-Brito et ak (2010[>, Johnson et al.| 


( |201 H >, Johnson et al. ( |2013| ), and Gonzalez et al. ( |2011) found 
a similarity between the [ar/Fe] abundance ratio of Bulge and 
thick disc giant stars. The origin of this discrepancy has been 
pointed to a systematic difference when analysing the abun¬ 
dances of dwarf stars from the disc and Bulge giants (Melendez 
et al. 2008 1 . In support of this scenario, the microlensed dwarfs 


from the Bulge have shown the same similarity in the a-element 
enhancement to those of the thick disc dwarfs (Bensby et al. 
201 1]>. Therefore, the usage of well-calibrated measurements on 


homogenous datasets are fundamental for these kinds of studies. 

A few high-resolution studies have looked at a few fields 
located at relatively similar latitude between each other (e.g. 
Johnson et al. 2013j>, but always at latitudes higher than b = -4°, 


where the contribution of the metal-rich stars is known to de¬ 
crease ( |Zoccali et al.||2008) |Johnson et al.|[20l3| |Ness et ak] 
2013a| l. The limited samples could therefore miss some subtle 
gradients. This is now investigated in this work with ~100 stars 
per field in a region where the two populations are well repre¬ 
sented. This study provides important constraints on the global 
view of the metallicity and cr-element distribution of the Bulge 
which are fundamental to reconstruct its formation scenario. 


2. GIBS high-resolution spectra 


The Giraffe Inner Bulge Survey (GIBS) is a survey of ~ 6,500 
RC stars in the Milky Way bulge observed with the GIRAFFE 
spectrograph of the FLAMES instrument (Pasquini et al .|2000 1 
at the ESO Very Large Telescope (VLT). The aim of the GIBS 
survey is to investigate the metallicity and radial velocity dis¬ 
tribution of Bulge stars across different fields, spread over a 
large area of the inner Galactic bulge (-10° < l < +10° and 
-10° < b < +4°). A detailed description of the target selection 
and data products of the program is given in Paper I. In this ar¬ 
ticle, we focus on the analysis of the high-resolution spectra of 

400 RC stars located in four fields at Galactic latitude b -3.5° 

for which we provide a brief description summary. 

Target stars for the spectroscopic observations were selected 


from the VVV multi-band photometric catalogues (Gonzalez 
et al.] |2012) >. The location of the fields was carefully selected 
in order to overlap with additional optical photometric observa¬ 
tions from the OGLEII survey (Sumi et al .|2004| l. This allows us 
to add extra information on the target stars such as proper mo- 




Fig. 1 . The location of the fields analysed in the present study 
(red circles) overplotted on an optical image of the Milky Way 
bulge (©Serge Brunier). The location of the fields discussed 
in previous studies (Zoccali et al. 2008 


Hill et al. 


2011 ) 


are 


shown as Z08 and Hll, as well as the location o f|Ness et al.| 
( |201 3a| >, [Johnson et al.| ( |20TT] i, and |Johnson et al.] ( |2013 1 fields 
as ARGOS, Jll and J13, respectively. Lower panel shows the 
colour-magnitude diagram for field HRm5m3. Red filled circles 
show the corresponding target stars for this field. 


tions and a larger colour baseline (V - K s ) for deriving the initial 
values for effective temperatures. 

Figure |T| shows the location of the 4 fields analysed in the 
present study overplotted on an optical image of the Milky Way 
bulge. The location of other fields with previous studies simi¬ 
lar to the one presented here (|Zoccali et al. 2008; Ness et al. 
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|2013a[ |Johnson et al.||20lTl |2013[ > are also shown. The colour- 
magnitude diagram for one of our observed fields (HRm5m3), 
where our selected targets are clearly highlighted, is also shown 
in Fig. |T] The selection box of targets in the colour-magnitude 
diagram can be clearly seen in Fig[I]and it was designed to tar¬ 
get RC stars while minimising the contamination for foreground 
disc stars. 

A description of the observations is reported in Table |T| All 
spectra in the fields analysed in this study were observed using 
the high-resolution grating (R=22,500) through setup HR 13 cen¬ 
tered at ~ 6300A in order to measure the chemical abundance of 
iron and a-elements in a similar way as in Gonzalez et al. ( 2011| . 

Reduction of the spectra was carried out using the ESO 
GIRAFFE pipeline. Bias and flat-held correction, individual 
spectral extraction, and wavelength calibration have been ap¬ 
plied by using the ESO/GIRAFFE pipeline^] An adequate sky 
spectrum for each exposure was obtained by co-addition of sky 
spectra obtained from dedicated fibres in each GIRAFFE con¬ 
figuration. Sky subtraction was carried out on each spectrum us¬ 
ing skytweak task in IRAF as described in Paper I. Heliocentric 
radial velocities were measured using IRAF/xcor task by cross¬ 
correlation against a synthetic spectra for a typical bulge K giant 
star (T eff = 4500K, log# = 2.3 and [Fe/H] = -0.3), cover¬ 
ing the corresponding wavelength range of the GIRAFFE HR 13 
setup from 6100A to 6400A. For each star, multiple exposures 
(~10) were individually corrected for radial velocity and indi¬ 
vidual products were mean combined to produce the final set of 
spectra for each held. The errors measured on radial velocities 
were typically ~ 1 km/s. 


3. Stellar parameters 

Stellar parameters, namely effective temperature (T e ff), surface 
gravity (logg), microturbulence velocity (£) and metallicity 
([Fe/H]) were derived based on an iterative method similar to 
the one described in Zoccali ~et al.| ( |2008| l and |Gonzalez et al. 
( 201 l| l. The method consists of the spectroscopic refinement of a 
set of the first guess photometric stellar surface parameters. The 
objective is to And the best combination of parameters in order 
to impose an excitation equilibrium of Fe lines (zero slope of Fe 
abundance as a function of excitation potential of Fe lines) and 
simultaneously measuring the same Fe abundance for all lines 
independently of their equivalent width (measuring a zero slope 
of Fe abundance as a function of log(EW//l)). Note that we rely 
on the usage of only Fel lines for this procedure, as the resolution 
GIRAFFE spectra does not provide us with sufficient clean Fell 
lines across our analysed spectral region. In Zoccali et al. ( 2008) 
and |Gonzalez et aT7| ( |201 1[ > this process was done manually, in 
the sense that the parameters were modified by hand to minimise 
the corresponding slopes. This manual procedure was found to 
produce some systematics, particularly in the high metallicity 
regime ( Hill et al.|201 1 1 . To correct for this issue and bring our 
datasets into a common baseline, we have now improved our 
method into an automatic procedure, similar to the one used in 
([Hill et al.||201l| ) but using the GALA code (Mucciarelli et al. 
2013) to spectroscopically refine the first guess photometric pa¬ 
rameters. Note that another difference with respect to the afore¬ 
mentioned studies is that for this work we have used ATLAS9 
model atmospheres (Castelli & Kurucz 2004). This is because 
GALA already includes a dynamic call to the ATLAS9 code so 
that whenever GALA needs to investigate a new set of atmo¬ 


spheric parameters, ATLAS9 is called and a new model atmo¬ 
sphere is computed. 

As a first step, the photometric temperature is calculated for 
each star using the (V - K s ) colours from OGLEII (Udalski 
2002] > and VVV survey catalogues (Gonzalez et al.|2012|7 


et al. 


dereddened based on the high resolution extinction maps from 


Gonzalez et al. (2012), and applying the Ramirez & Melendez 


(2005) calibration. Absolute V band magnitudes, calculated 


from the measured distances to each specific field from Gonzalez 
|et al. (j2013| and bolometric corrections from Alonso et al. 


(1999 1 , are then used to estimate photometric gravities based on 


1 http://www.eso.org/sci/software/pipelines/ 


the usual formula: 

log (g) = log (go) + log (^4) + 0.4 (M BoU - M Bo i,o) + 4 log 

where M Bo i,o = 4.72, T e fr, 0 = 5770 K, and log(g Q ) = 4.44 
dex. A fixed value of M, = 0.8M o has been adopted similarly 
to |Zoccali et al.| ( |2008| ) and |Hill et al.| ( [20TT] ). Microturbulence 
velocity and global metallicity are set to 1.5 and 0.0, respec¬ 
tively, as a first step. These values are used to obtain a first 
guess ATLAS9 stellar model atmosphere and are subsequently 
refined spectroscopically by GALA using the equivalent widths 
of isolated Fe lines obtained by means of DAOSPEC (Stetson & 
|Pancino|2008| ). 

After feeding GALA with the photometric stellar surface pa¬ 
rameters and the corresponding model atmosphere, the code iter¬ 
atively searches for spectroscopic effective temperatures and mi¬ 
croturbulence velocity by imposing excitation equilibrium and 
the null slope of iron abundance versus equivalent width of the 
Fe lines. During each iteration, a new model atmosphere is gen¬ 
erated by GALA using the refined stellar parameters. Although 
available in GALA, we do not refine log g values by requiring 
ionisation equilibrium, but we rely on the values derived photo¬ 
metrically. As discussed in Zoccali et al.| (2008l, the resolution 
of GIRAFFE is not high enough to resolve a sufficient number 
of clean Fell lines in our spectra and thus using the few avail¬ 
able lines would introduce more errors than just using the pho¬ 
tometrically derived values for logg. Therefore, after a first set 
of best-fitting values for T c n , /■’, and [Fe/H] is found by GALA, 
we re-calculate the logg value according to the new parameters 
and perform another set of GALA iterations to refine the stellar 
parameters with the new log g. 

The list of Fel lines used for this method is the one used in 
Lecureur et al. (2007) following a careful cleaning of blended 
lines. The check for blends of the Fe lines was done using 
a spectrum of /r Leo obtained at the Canada-France-Hawaii 
Telescope with the ESPaDOnS spectro-polarimeter at a resolu¬ 
tion R -80,000 and high S/N per pixel (~500). As a first step, the 
resolution of the /j Leo spectrum was degraded to the resolution 
of GIRAFFE HR 13 setup of R-22,500. This was followed by a 
line-by-line inspection performed over the wavelength range of 
our science data. Several blended Fe lines were removed from 
the list ending with a total of 27 non-blended Fe lines in the 
HR13 setup. The log gf values from the original line list from 
Lecureur et al.|< j2007) > are calibrated so that each line provides an 
abundance of 0.30 dex from the EW measured in the spectrum of 
/j. Leo while adopting T e ff = 4540 K, log g = 2.30, and /' = 1.30. 
Therefore, in order to further check our final Fe line list, we mea¬ 
sured the EW of the corresponding line using DAOSPEC in the /j 
Leo spectrum (now at GIRAFFE resolution) and used GALA to 
derive the Fe abundance by fixing the model to those same stel¬ 
lar parameters from Lecureur et al. ( 2007| ). GALA retrieved an 
abundance of [Fe/H] =0.32, in excellent agreement with the lit¬ 
erature value of [Fe/H]=0.34 ( [Gratton & Sneden|1990| ). On the 
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-1.5 -1.0 -0.5 0.0 0.5 

[Fe/H] (this work) 


Fig. 2. Difference between metallicity (upper panel), effective 
temperature (middle panel) and microturbulence velocity (lower 
panel) for a sample of RC (red tilled circles) and RGB (black 


empty circles) derived in this work and as derived in Hill et al. 
( |2011 [ >. Solid lines in each panel indicate the mean difference of 
each sample and the dashed lines shows the lcr range around the 


other hand, by leaving the stellar parameters free with the excep¬ 
tion of log g which was fixed to log g = 2.30, GALA retrieves 
T eff = 4490 K, £ = 1.43 and [Fe/H 1=0.31 for /j Leo. These val¬ 
ues are in good agreement with those reported in the literature 
( jSmith & Ruck|2000} |Gratton & Sneden|1990) >. The final list of 
Fel lines used in our analysis is provided in Table [2] 

Particular care has been devoted to ensure that the adopted 
stellar parameters were i) sufficiently accurate to allow reliable 
abundance estimates, and ii) consistent with previous analysis to 
guarantee a homogeneous comparison across various samples. 
This provides us with a test for our stellar parameter measure¬ 
ments while at the same time ensures that our newly derived 
abundances can be safely compared to our previous work. With 
this aim, we have re-derived the abundances for the Baade’s 


Window sample of Red Giant Branch (RGB) stars from Zoccali 
et al. ( j2008| ) and for the RC sample from Hill et al. ([201 I) using 


the GALA-based procedure described above. Figure[2]shows the 
comparison between the stellar surface parameters and [Fe/H] 
for RC and RGB samples from Hill et al. ( 2011| and those de¬ 
rived with our method. This comparison involves only the calcu¬ 
lation of stellar parameters, since the EW measurement is done 
using DAOSPEC both in |Hill et ak ( |201 1[ ) and in this work. 

As shown in Fig. [2] there is a good agreement between both 
methods with no significant differences for either of the two 
samples analysed (RGB and RC stars). In particular, for effec¬ 
tive temperature we found a mean difference AT e ff = 2 + 199 K 
and -23 + 185 K for RGB and RC, respectively. On the other 
hand, micro-turbulence velocities result in a mean difference of 
A£ = -0.14 + 0.19 km/s for RGB stars and A£ = -0.12 + 0.26 
km/s for RC stars. These stellar parameters translate into a mean 
metallicity difference of [Fe/H] = -0.02 ± 0.19 and [Fe/H] = 
-0.01 + 0.15 dex. The derived differences are negligable and the 
scatter between both samples is in good agreement with the ex¬ 
pected errors on the measurements of the surface stellar parame¬ 
ters as discussed in|Hill et ak|(j2011|. Furthermore, this confirms 



-1.5 -1.0 -0.5 0.0 0.5 

[Fe/H] (this work) 


Fig. 3. Difference between the [Fe/H] abundances obtained with 
the adopted GALA procedure using the correct set of initial 
stellar surface parameters and using fixed values (T e ff=4500 K, 
£=1.5 km/s, log g= 1.9, and [Fe/H | =0.00) for the RC sample 
from [Hill et ;iT ( 20T~I ) 


that our results and those from Hill et al. 

(2011) are in good 

agreement, with no systematic differences 

tetween metalicities 

obtained for RGB stars from Zoccali et al. ([2008) and for RC 

stars as previously reported in 

Hill et al. 

(2011). Therefore, the 


measurements on the Baades Window samples (RC and RGB) 
can be compared in a consistent and homogenous way with the 
results from the GIBS survey discussed in the present study. 

Finally, we investigated the dependence of the final [Fe/H] 
abundance on the initial photometric estimations of stellar pa¬ 
rameters. Reddening corrections and distance spread trans¬ 
late into uncertainties in the initial photometric parameters. 
Furthermore, adopting an effective temperature calibration dif¬ 
ferent to the one of [Ramirez et al.| ( |20001 > would cause a differ¬ 
ent starting point for the first run of GALA. In order to evalu¬ 
ate the impact of these parameters in the final abundances, we 
have tested a run of the RC sample from Hill et al. ( 2011) >, 
but starting from the same set of stellar parameters for all the 
stars: T e ff=4500K, £=1.5 km/s, logg=1.9 dex, and [Fe/H]=0.00. 
Figure [3] shows the comparison between the original [Fe/H] 
for this sample and the one resulting from the same stellar pa¬ 
rameters. Because the mean difference in the [Fe/H] estimates 
is A[Fe/H]=0.02 dex, no significant systematic offsets are ex¬ 
pected in the final set of [Fe/H] abundances due to the adopted 
set of initial stellar parameters, no matter how far they are from 
the correct values. On the other hand, the dispersion between the 
[Fe/H] is found to be o-[ Fe / H ] =0.13 dex, which is not negligible. 
Differences of up to 0.3 dex can be found among the most metal- 
rich stars as a result of using an incorrect set of initial parame¬ 
ters. It is worth mentioning that GALA has an optional initial 
stage, called the guess working block , where the stellar parame¬ 
ter space is explored in a coarse grid in order to verify and refine 
poorly known initial parameters. However, in our case concern¬ 
ing Bulge giant stars, the reddening and distance spread uncer¬ 
tainties are always propagated into the calculation of the photo¬ 
metric surface gravity. The non-variation of this value against the 
procedure of GALA makes the guess working block redundant 
and an error similar to the A[Fe/H]=0.13 dex observed in this 
test is expected within the uncertainty of this particular parame¬ 
ter. Most likely, the observed uncertainty associated to the initial 
stellar parameters is not caused by GALA itself, but for the in¬ 
ability to optimise the entire grid of parameters simultaneously 
in the bulge. Were a larger number of Fel and Fe II lines avail- 
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able in our spectra, the final scatter on stellar parameter would 
converge to less than 50 K in T e g- and 0.09 in log g as described 
in (Mucciarelli et al7] |2013[ ) thus decreasing the relatively large 
scatter in [Fe/H] observed in Fig. [3] 

The uncertainties in the stellar parameters are derived inter¬ 


nally by GALA. As pointed out in Mucciarelli et al. ( 2013[ ), the 
errors in the stellar parameters are dominated by the number of 
Fel lines used and their distribution in both ^f ex range and tran¬ 
sition values. For each stellar parameter, an optimisation param¬ 
eter is calculated by applying a jackknife boot-strapping tech¬ 
nique. This means that if the optimisation is computed using a set 
of N spectral lines, the parameters will then be re-calculated N 
times, but using N-l lines each time. This calculation is then re¬ 
lated to a parametrisation factor that measures the way each pa¬ 
rameter affects the slopes and iron abundance differences within 
the GALA procedure. This parametrisation value is calculated 
by GALA by varying the best value found around the local min¬ 
imum. This method and their respective equations are described 
in detail in Mucciarelli et al. ( 2013| . In our sample, the resulting 
errors are found to be <XT dI =245 +99 K, cr^- 0.29 ±0.11 km/s. 
Each parameter is then varied by its corresponding uncertainty, 
while keeping the others unchanged, and the abundances are re¬ 
derived. The individual effects in abundance from the variations 
of each parameter are added in quadrature to obtain the final 
error in metallicity. In our sample, this computation results in 
a mean metallicity error of cr[ Fe / H ]=0.20 ± 0.07 dex. However, 
note that the errors in [Fe/H] can reach up to 0.4 dex for the 
most metal-rich stars but can be lower than 0.1 dex for metal- 
poor stars and are mostly driven by T e ff uncertainties. These er¬ 
rors are consistent with those reported in |Hill et~ak ( 201 1[ > and 
Zoccali et al. ( |2008| > for similar datasets and are also in agree¬ 
ment with the differences seen in Fig. [2] Furthermore, the larger 
errors in stellar parameters towards the high metallicity regime 
are expected. Extensive analysis of errors in derivation of stellar 
parameters was presented in Smiljanic et al. (2014) in the con¬ 
text of the Gaia-ESO survey and showed how the errors become 
larger at high metallicities due to the increase of line blends. 


4. Alpha element abundances 

In order to measure a-element abundances for all the target stars 


we adopt the same procedure and line lists used in Gonzalez 
et al. (2011|). As a first step, the corresponding ATLAS9 model 


atmosphere (Castelli & Kurucz 2004) is generated using the stel¬ 
lar parameters found by GALA for each star. This model is then 
fed to MOOG (v. Feb2013, Sneden|1973| l, in order to generate 
a synthetic spectrum. The synthetic spectrum is iteratively com¬ 
pared with the observed one after a local normalisation of the 
continuum, by varying only the elemental abundance of interest 
in each iteration until the best fitting abundance if found. 

Based on similar dataset - same setup, resolution and simi¬ 
lar signal-to-noise - Gonzalez et al. ( 2011[ ) showed that among 
all the o'-clements measurable in the observed spectral range Mg 
has the smallest scatter. Therefore, we consider Mg to be the key 
element allowing a precise comparison between the abundance 
trend observed in different bulge and disc star samples. However, 
in order to obtain precise abundances of Mg we need to first cal¬ 
culate Ca abundances. Indeed, as discussed in detail i n|Lecureur| 
|et al,| ( [2007| >, |Hill et al.| ( |201 1[ >, and |Gonzalez et al7| ( |20 11 ] > the 
continuum around the Mg triplet at 6319A is severely affected 
by 6318.1A Ca I autoionisation line. The large broadening of 
this Ca line causes a drop in the Mg triplet region that can result 
in an incorrect abundance measurement for Mg. For this reason, 
we first measured Ca abundances using two lines at 6162.1 and 


0GLE26-415574: [Fe/H]=0.15, [Co/Fe]=0.09, [Mg/Fe]=0.06 



6316 6317 6318 6319 6320 6321 6322 

Wavelength 


Fig. 4. Observed (black circles) and synthetic spectra for a star 
(OGLE26-415574) in the field HRm5m3 zoomed into the region 
used to derive Ca (upper panel) and Mg (lower panel) abun¬ 
dances. The synthetic spectra created with the best fitting Ca 
([Ca/Fe]=0.06) and Mg ([Mg/Fe]=0.05) abundances are shown 
in solid red lines. The elements producing the stronger absorp¬ 
tion lines are marked in each panel. For each spectral region we 
show a zoomed-in section around the Ca line at 6166.4 A and the 
Mg triplet that includes the synthesis produced using the best- 
fitting abundance as red solid lines. In blue and green we show 
the synthetic spectra using a 0.1 dex abundance variation around 
the best-fitting value. The Y-axis of both zoomed-in regions have 
the same scale with normalized fluxes that range from 0.8-0.5 
for Ca and 1.0-0.7 for Mg lines 


6166.43A (see Fig. [f] lower panel and Table [2j. Once the Ca 
abundances have been measured we provide this value as an ad¬ 
ditional input for MOOG to be included in the synthesis used 
for the calculation of Mg abundances. The local normalisation 
of the continuum is then manually improved as an additional 
way to account for any failure to reproduce the underlying Ca 
autoionisation line. An example of this fit is shown in Fig. [4] 
This procedure was applied to all the analysed stars, obtaining 
Ca and Mg abundances for a total of 400 stars spread in our four 
fields. Fig.[4]also shows an example of the sensitivity of the lines 
to the variations in [Mg/H] and [Ca/H]. For a variation of +0.1 
dex, the effect in the fit region of Mg is clearly seen, while the 
Ca lines appear to be less sensitive to such variations. Both of 
the Ca lines that are used in this work show a similar behaviour. 
The difference in [Ca/Fe] abundances measured from each of 
the lines has a cr of 0.05 dex, which is negligible compared to 
the errors involved in our analysis due to the stellar parameters 
uncertainties. 

The uncertainties on the measured [Mg/H] and [Ca/H] abun¬ 
dances were calculated in the same way as in Gonzalez et al. 
( 201 1| > by varying the stellar parameters by their correspond¬ 
ing uncertainties and re-calculating the best fitting abundances. 
These errors are then added in quadrature together with the error 
from our spectral synthesis fitting procedure of 0.1 dex derived 
i Gonzalez et al. (201 1)>. The errors are individually reported in 


Tableland are found to be of the order of 0.14 dex for [Mg/H] 
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Table 2. Atomic line list for Mg and Ca using in this work. Also 
listed are the excitation potential (Tex) and oscillator strength 
(log gf) for each analysed line. 


A (A) 

Element 

log gf 

/Tex 

6120.246 

Fel 

-5.970 

0.915 

6137.691 

Fel 

-1.375 

2.588 

6151.617 

Fel 

-3.312 

2.176 

6157.728 

Fel 

-1.160 

4.076 

6162.160 

Cal 

-2.720 

1.899 

6165.360 

Fel 

-1.470 

4.143 

6166.430 

Cal 

1.142 

2.521 

6173.334 

Fel 

-2.880 

2.223 

6180.203 

Fel 

-2.617 

2.727 

6187.989 

Fel 

-1.620 

3.943 

6191.558 

Fel 

-1.416 

2.433 

6200.312 

Fel 

-2.405 

2.608 

6213.430 

Fel 

-2.481 

2.223 

6219.281 

Fel 

-2.434 

2.198 

6226.734 

Fel 

-2.120 

3.883 

6229.226 

Fel 

-2.805 

2.845 

6230.722 

Fel 

-1.279 

2.559 

6246.318 

Fel 

-0.805 

3.602 

6252.555 

Fel 

-1.727 

2.404 

6253.829 

Fel 

-1.299 

4.733 

6271.278 

Fel 

-2.703 

3.332 

6301.500 

Fel 

-0.718 

3.654 

6311.500 

Fel 

-3.141 

2.831 

6318.710 

Mgl 

-2.000 

5.110 

6319.230 

Mgl 

-2.240 

5.110 

6319.490 

Mgl 

-2.680 

5.110 

6322.685 

Fel 

-2.448 

2.588 

6330.848 

Fel 

-1.64 

4.733 

6335.330 

Fel 

-2.177 

2.198 

6336.823 

Fel 

-0.856 

3.686 

6355.028 

Fel 

-2.32 

2.845 

6380.743 

Fel 

-1.475 

4.186 


and 0.21 dex for [Ca/H]. Figure [5] shows the error estimation for 
Ca and Mg abundances as a function of the estimated uncertain¬ 
ties for the stellar parameters of each star. It can be clearly seen 
that Ca abundances are much more sensitive to any variations on 
stellar parameters than Mg abundances. Furthermore, the larger 
sensitivity of Ca abundances to stellar parameters is seen at all 
metallicities, but is even larger in the metal-rich regime. 


5. Results 

Stellar parameters, [Fe/H], [Ca/Fe] and [Mg/Fe] values mea¬ 
sured for a total of 400 stars in the GIBS sample are listed in 
Tableland discussed in detail in this section. 

5.1. The Bulge metallicity distributions at constant latitude 
The derived metallicity distribution of the five fields located at 



0 100 200 300 400 500 

” Cell) 
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[Fe/H] 


Fig. 5. Abundance measurement errors as a function of errors in 
the stellar parameters: [Fe/H] (upper left panel), microturbulence 
velocity (upper right panel) and effective temperature (lower left 
panel). Errors as a function of the corresponding [Fe/H] value 
are also displayed (lower right panel). In all panels errors for 
[Ca/H] abundances are shown as red empty circles and errors 
for [Mg/H] as red filled circles. 


et al. 1998). We have adopted a Gaussian kernel and the optimal 
smoothing parameter as defined in Silverman et al. (1998). 

The good number statistics of the total sample allows us to 
investigate in detail the metallicity distribution. In particular, we 
are interested in evaluating the decomposition of the distribution 
in two populations, as suggested in previous studies (|Hill et al. 


201 l[|Roja s-An~iag ada et al.|2014 1 . Note that Ness et al.| ( 20l3a| l 


suggested that the Bulge metallicity distribution can be decom¬ 
posed into five Gaussian components, each of them correspond¬ 
ing to different populations from the bulge and foreground disc. 
Here, we limit this decomposition to only two Gaussian compo¬ 
nents. The main reason for setting this limit is because studies, 
based on both bulge kinematics and morphology, have found evi¬ 
dence for only two main components: a metal-poor spheroid-like 
component and a metal-rich bar-like component. Also, in this 
study we are restricted to inner Bulge fields in which we have 
a narrow target selection box in the colour-magnitude diagram. 
For this reason, the contamination by the thin/thick disc and halo 
is negligible here (components D, F, and E in Ness et al.|2013a) >. 
It is worth stressing that here the use of two Gaussian distribu¬ 
tions is not necessarily meant to describe the actual shape of the 
metallicity distribution of each of the two components. It repre¬ 
sents the attempt to perform the best possible parametrisation of 
the observed global metallicity distribution, so we can therefore 
investigate how its shape changes across the fields. 

We perform a least-squares Gaussian fit to the global metal- 



statistics, we initially construct the global metallicity distribu- 


ifying an initial guess for the mean, peak, and sigma value. 


tion from all the stars in the 5 fields together. Fig. [6] shows the 


The code then refines the fit for both Gaussian functions si- 


resulting global metallicity distribution. In order to avoid bin¬ 
ning effects when looking at the shape of the distribution, we 
have estimated the probability density underlying in the global 
sample by using the kernel density estimator method (Silverman 


multaneously. The best-fit Gaussian distributions and the global 
metallicity distribution are shown in Fig. [6] The resulting pa¬ 
rameters for the best-fit metal-poor Gaussian component are 
< [Fe/H] >=-0.31 and cr[Fe/H]=+0.31. On the other hand, the 


7 




































O. A. Gonzalez et al.: The Giraffe Inner Bulge Survey (GIBS) II 


Table 3. Stellar parameters, [Fe/H], and [Mg/Fe], [Ca/Fe] abundance ratios for stars in the four high-resolution GIBS fields. The 
full table is available in electronic form. Individual coordinates and the associated photometry is available in Zoccali et al. (20141. 


Field 

STAR 

T e ff 

0" T eff 



logg 

[Fe/H] 

cr[Fe/H] 

[Ca/Fe] 

cr[Ca/Fe] 

[Mg/Fe] 

°"[Mg/Fe] 

HRm5m3 

OGLE26_649743 

4328 

240 

1.6 

0.30 

2.29 

0.14 

0.24 

0.02 

0.21 

0.14 

0.09 

HRm5m3 

OGLE26_661709 

4537 

269 

1.4 

0.25 

2.22 

0.28 

0.33 

0.01 

0.19 

-0.04 

0.09 

HRm5m3 

OGLE26_649746 

4501 

257 

1.0 

0.24 

2.36 

0.35 

0.21 

0.15 

0.20 

0.01 

0.16 

HRm5m3 

OGLE26_661748 

4468 

151 

1.3 

0.28 

2.52 

0.24 

0.19 

0.00 

0.22 

0.01 

0.21 




[Fe/H] 


Fig. 6. Overall metallicity distribution obtained from the com¬ 
bination of the 4 GIBS high-resolution fields and the RC sam¬ 
ple of Hill et al. (2011). The probability density distribution is 
overplotted in the upper panel as a dashed line. The correspond¬ 
ing variance on the calculated probability densities is shown in 
blue. The lower panel shows the best two Gaussian fit to the up¬ 
per panel distribution. The resulting distribution from the two 
Gaussians is shown as a solid line. 


metal-rich Gaussian component has a mean < [Fe/H] >=+0.26 
and cr[Fe/H]=+0.20. The two Gaussian fits to our global metal¬ 
licity distribution are in good agreement with those of FfiTi] 
et al. < 2011] ) who finds a metal-poor component with mean 
of [Fe/H]=-0.27 and a metal-rich one of [Fe/H]=+0.29. The 
sigma of the metal-poor component is in good agreement with 
that of |Hill et ah] ( j201 1) >, finding a broad metal-poor compo¬ 


nent with <x[Fe/H]=0.31. On the other hand, [HTlI et al.||2011j l 
find a metal-rich component slightly narrower than ours with 
cr[Fe/H]=0.12. Nevertheless, the results are fully compatible. 
The observed difference in the width of the metal-rich compo¬ 
nents originates from the increase of the [Fe/H] errors towards 
the metal-rich regime and the respective error deconvolution that 
was applied in the Hill et al. poll} analysis. 

The individual metallicity distributions for all the fields pre¬ 
sented in this work are shown in Fig. [7] From the measured 
mean metallicity in each field we see no evidence of a gradi¬ 


ent such as the one seen along the minor axis ( Zoccali et al. 
2008; [Ness et aT7||2013a( >. The observed mean [Fe/H] values of 
the fields studied in this work range between [Fe/H] =0.10-0.20, 
consistent with what is shown in the mean photometric metal¬ 
licity map from Gonzalez et al. ( [2013) . We do not attempt a 
multi-Gaussian fit to the individual field metallicity distributions 
because we do not have enough statistics in each field, there¬ 
fore preventing us from deriving a statistically robust conclu¬ 
sion. However, we use the probability density distribution that 
was obtained from the global metallicity distribution as a refer¬ 
ence in order to evaluate field-to-field variations. From Fig. [7] we 
notice that the metallicity distribution of individual fields span 
the same range: — 1 < [Fe/H] < +1. The general shape of the 
distributions is well represented by the probability density dis¬ 
tribution constructed over the global metallicity distribution. A 
more noticeable difference can be observed in field m5m3 where 
the metal-rich side of the distribution shows a strong peak near 
Solar metallicity that is not seen in the other fields. 

In order to evaluate whether the metallicity distributions for 
each field originate from different populations or not, we per¬ 
formed a two sample Kolmogorov-Smirnov (K-S) test and com¬ 
pared the metallicity distribution found in Baade’s window with 
the individual GIBS field distributions. The K-S test returns p- 
values (i.e. the probability that two functions belong to the same 
distribution): 0.22 for p8m3, 0.48 for p4m3, 0.03 for m5m3, 
and 0.16 for m7m3. With the exception of field m5m3, the K- 
S test confirms that our samples are most likely drawn from the 
same parent population. For field m5m3, the result is not con¬ 
clusive. From the visual inspection of the metallicity distribution 
we notice that the main difference with respect to the other fields 
comes from a concentration of stars at [Fe/H]—(-0.1. 

We conclude in favour of a null mean metallicity gradient 
along the Bulge major axis, as previously suggested in Johnson 


et al. (201 3)> at higher latitudes (but see Rangwala & Williams 


( |2009] > and |Babusiaux et ak| ( |2014| l). Furthermore, we conclude 
that the metallicity distributions presented in this work for fields 
along the longitude strip at b=-4° originate from the same 
parent population. Based on this, we suggest that the relative 
contribution of the potentially two different Bulge components 
(spheroid-like and bar) does not change along the major axis, 
providing that this is indeed the origin of the metallicity gradient 
found along the minor axis. 
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Fig. 7. Metallicity distributions for the four GIBS high-resolution fields. Additionally, the metallicity distribution for Baade’s win¬ 
dow from Hill et al. (2011) is shown in red in the middle panel. For each of the metallicity distributions we overplot the probability 
density distribution for the entire sample of metallicities shown in Fig. [6] 


5.2. The Bulge alpha-element abundances at constant 
latitude 


We now present the results for [Ca/Fe] and [Mg/Fe] abundances 
for all our target stars, obtained by comparison to synthetic spec¬ 
tra created using the corresponding stellar surface parameters 
and [Fe/H] discussed in previous sections. A total sample of 
400 stars have been analysed and their [Ca/Fe] and [Mg/Fe] over 
[Fe/H] trends are shown in Fig. [8] Also in Fig.[8]we show a com¬ 
parison of our results with respect to those derived from a similar 


sample in Baade’s window measured in Gonzalez et al. (2011!. 


While the scatter of the measurements in [Ca/Fe] are com¬ 
parable in both samples, the [Mg/Fe] trend shows a remarkable 
reduction of the scatter with respect to the previous measurement 
from Gonzalez et al. ( |2011| >. This is most likely a consequence 
of our improved calculation of stellar parameters as discussed in 
Section 3. The larger dispersion in Ca abundances with respect 
to Mg can be clearly seen in Fig. [8] This spread in [Ca/Fe] is 
most likely caused by the larger dependence of Ca abundance 
on stellar parameter uncertainties. In Fig. [9] we show that when 
only the stars with small uncertainties are considered, the se¬ 
quence of [Ca/Fe] becomes tighter, with much less scatter at all 
metallicities. It can be seen that stars with high [Ca/Fe] found 
at high metallicities [Fe/H] > 0 are restricted to those measure¬ 
ments with large uncertainties. In our measurements any error 
in Ca abundances should in principle translate in an error in Mg 
abundances due to the underlying auto-ionisation Ca line that 
strongly affects the continuum around the Mg triplet (Hill et al. 
2011) >. Therefore, the adoption of an incorrect [Ca/Fe] abun¬ 
dance would result in a larger scatter for [Mg/Fe] abundances 
due to a poorer quality fit in the respective region. We do not see 
this large spread in Mg abundances, most likely because this is 
partially solved by the local normalisation of the continuum. As 
a matter of fact, the resulting [Mg/Fe] trend for the Bulge can be 
traced with a great level of detail. 

We compare the trends of [Ca/Fe] and [Mg/Fe] as a function 
of [Fe/H] for our entire sample with those from giant stars of 
the local thin and thick disc published in Gonzalez et al. ( |2011) , 
As seen in both Ca and Mg abundances in Fig. [8] the a-element 
enhancement of the Bulge with respect to the thin disc is found 
at all metallicities. The thick disc on the other hand appears as 
enhanced in a-element abundances as the metal-poor Bulge gi¬ 
ants, though covering a much smaller range in [Fe/H]. These 
patterns in the a-element abundances of Bulge stars with re¬ 
spect to the disc abundances has been already observed in sev¬ 
eral studies, based on different samples and measurement tech 


niques (Melendez et al. 2008 [ 

Alves-Brito et al. 2010 

Gonzalez 

et al.EOl 1; Bensby etal.2010 

2011 Ryde etal.2010; 

Hill et al. 



[Fe/H] 


Fig. 8. Distribution of [Mg/Fe] (upper panel) and [Ca/Fe] (lower 
panel) as a function of [Fe/H] used as a diagnostic of the for¬ 
mation time-scale of the Bulge. Shown as red circles are the 
[Mg/Fe] and [Ca/Fe] values obtained from the HR 13 spectra of 
400 stars in the four GIBS fields presented in this work com¬ 
pared to those for RGB stars from the Baade’s window sample of 
Gonzalez et al. (2011) which are shown as black empty circles. 
Measurements for the local thin (black squares) and thick (blue 
triangles) disc stars from the sample from Gonzal ez et al.| < 2011 1 
are also shown in both panels. The uncertainties for [Ca/Fe] 
and [Mg/Fe] abundances are shown as a light-blue contour in 
both panels according to the individual uncertainties reported in 
Table [3 


|20 11 [ [Johnson etal.|2011||2013||2014| l. In addition, we have cal¬ 
culated the metallicity at which the [Mg/Fe] abundances starts 
to decrease (the so-called knee) similarly as in Gonzalez et al. 
( 2011| >, by producing a linear fit to the metal-poor ([Fe/H]=-1.6 
to -0.6) and metal-rich ([Fe/H |=-0.3 to +0.3) stars. The inter¬ 
section between both linear fits corresponds to the a-element 
knee located at [Fe/H] =-0.44, in good agreement with the val¬ 
ues found in |Gonzalez et al.| ( |2011] ) in a field located at the same 
latitude b=-4°. 


9 























































O. A. Gonzalez et al.: The Giraffe Inner Bulge Survey (GIBS) II 



0.4 


0.2 - 


0.0 - 


- 0.2 - 



• Bulge GIBS 
A Bensby 2014 
□ Jonhson 2014 


0.4 - 


TT 0.2 - 


0.0 


- 0.2 



- 1.5 - 1.0 - 0.5 0.0 

[Fe/H] 


0.5 


Fig. 9. Distribution of [Mg/Fe] (upper panel) and [Ca/Fe] (lower 
panel) as a function of [Fe/H] as in Fig. 11 Here we colour- 
coded the abundances of [Ca/Fe] and [Mg/Fe] according to the 
estimated uncertainty value. The colour-coding was fixed in both 
panels to the range of [Mg/Fe] uncertainties between 0.07 - 0.20 
dex from black to red. The uncertainty contour of [Mg/Fe] is also 
shown in the background of both panels as a reference. 


Bensby et al. ( |2013[ ) suggested that the position of the knee 
in [a/Fe] plot is found at 0.1-0.2 dex higher metallicities in the 
Bulge than in the thick disc. A similar result was proposed by 
Johnson et al. ( j2014f who added the analysis of Fe-peak ele¬ 
ments finding in particular that Co, Ni, and Cu appear enhanced 
compared to the disc. It is important to recall that the results 
presented in Johns orTet al.| ( 20 14) have been obtained by com¬ 
paring Bulge giants to dwarf stars from the disc. This technique 
might suffer from systematic offsets (Melendez et al. |2008 1 and 
it thus needs to be further investigated. The results from Bensby 
et al. ( 2013} do not suffer from these systematics, as the com¬ 


parison involves only dwarf and sub-giant stars for both disc 
and bulge; however, because the microlensed dwarfs sample is 
still relatively small, for the time being the derived result is far 
from being statistically robust. Unfortunately, in our work we do 
not have measured abundances for a sufficiently large number of 
thick disc giant stars to calculate the position of the correspond¬ 
ing [a/Fe] knee with enough accuracy to compare it with the 
value ([Fe/H] =-0.44 dex) measured in the GIBS fields. However 
we note that for the Bulge we find a [a/Fe] knee ~0.1 dex more 
metal-poor than the value presented in Bensby et al. (2013). 
This means that, assuming that our study is in the same scale as 
Bensby et al.( 2013[ ), our abundances would not favour a more 


metal-rich knee for the Bulge than the local thick disc. 

In Fig. [TO] we show a comparison between the [Mg/Fe] 
and [Ca/Fe] trends between GIBS giants, the giants analysed in 
Johnson et al/]( 2014] ), and the microlensed dwarfs from Bensby 


et al.|(|2013|l. The trends are very similar for [Mg/Fe] abundances 


while for [Ca/Fe] there is a clear systematic offset of the 
abundances with respect to both GIBS and . 


et al. 

(2013 

et al. 

(2014 


Bensby 


ohnson 


(2014). The fact that the offset is not seen in [Mg/Fe] indi- 


Fig. 10. Distribution of [Mg/Fe] (upper panel) and [Ca/Fe] 
(lower panel) as a function of [Fe/H] for the GIBS high- 
resolution fields studied in this work and for the measurements 
from Johnson et al. (2014) for giant stars (black squares) and mi¬ 
crolensed dwarfs from Bensby et al. ( 2013) ) (blue triangles). In 
the upper panel, the solid lines show the best-fit relation between 
[Mg/Fe] and [Fe/H] for the metal-poor and metal-rich regime. 
The intersection point is shown and represents the position of 
the knee. The uncertainties for [Ca/Fe] and [Mg/Fe] abundances 
of this work are shown as a light-blue contour in both panels. 


cates that most likely it is not due to a different method of spec¬ 
tral analysis, but it emphasizes the risk of comparing abundances 
for stars in very different temperature and/or gravity regimes (gi¬ 
ants versus dwarfs, in this case) because the abundances of some 
elements may be affected by different systematic errors. 

Finally, in Figure m we show the individual trends of 
[Mg/Fe] and [Ca/Fe] abundances as a function of [Fe/H] for 
the four fields presented in this work, compared to the [Mg/Fe] 
abundances for the local thick and thin disc. No evident field- 
to-field variations can be seen and the same similarity between 
the metal-poor Bulge (i.e. [Fe/H]< -0.5) and the local thick disc 
is found in all fields. This result for a strip at constant Galactic 
latitude is in agreement with the conclusions of Gonzalez et al. 
( |201 1| ) for the minor axis. 

6. Conclusions 

In this work we presented the surface stellar parameters, and the 
[Fe/H], [Ca/Fe], and [Mg/Fe] element ratios for 400 RC stars 
in four fields of the GIBS survey, corresponding to the high- 
resolution sample at b = -4° observed with GIRAFFE HR 13 
setup. Spectroscopic stellar surface parameters were obtained 
with an automatic method based on the usage of the GALA 
code. These parameters were used to obtain [Ca/Fe] and [Mg/Fe] 
abundances for each star based on fitting of synthetic spectra. 
The global and individual field metallicity distributions were 
constructed, in order to investigate the presence of gradients and 
field-to-field changes in the shape of the distributions. We also 
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[Fe/H] 


Fig. 11. Distribution of [Mg/Fe] (upper panel) and [Ca/Fe] 
(lower panel) as a function of [Fe/H] using different symbols 
for the four GIBS high-resolution fields studied in this work. 
The solid lines show the best-fit relation between [Mg/Fe] and 
[Fe/H] for the metal-poor and metal-rich regime obtained us¬ 
ing the entire sample. The uncertainties for [Ca/Fe] and [Mg/Fe] 
abundances are shown as a light-blue contour in both panels. 


investigated the [Ca/Fe] and [Mg/Fe] abundances as a function 
of [Fe/H] as a constraint for the time-scale of formation for the 
Bulge. 

The results can be summarised as follows: 

- A null metallicity gradient was found along Galactic longi¬ 
tude for a constant Galactic latitude. The metallicity distribu¬ 
tions of all the fields presented here are consistent with being 
originated from the same parent population. 

- A two Gaussian component provides a good fit to the global 
metallicity distribution at b=-4° with a narrow metal-rich 
([Fe/H]—t-0.26) population of stars and a broad metal-poor 
([Fe/H]—0.31) component. This is in agreement with previ¬ 
ous studies suggesting a bimodal population for the Bulge. 
This exercise has been carried on because several indepen¬ 
dent observations, including stellar kinematics and distance 
distribution, suggest the presence of two stellar populations 
in the Galactic bulge. Whether the two components have in¬ 
deed a different formation history, needs to be seen by model 
comparisons to observations as the ones presented here. 

- A very tight relation was found between [Mg/Fe] and [Fe/H], 
showing enhancement with respect to the solar ratio in the 
metal-poor regime, similar to the one seen in the local thick 
disc, and a decrease in [Mg/Fe] starting at [Fe/H]—0.44. All 
studied fields share these properties, with no appreciable dif¬ 
ferences among each other. The observed position of the so- 
called knee is in agreement with our previous results but it 
occurs at ~0.1 dex lower metallicity than what is observed in 
the bulge microlensed dwarf and sub-giants. 

- The [Mg/Fe] and [Ca/Fe] results presented here are fully 
compatible with previous results based on giant stars; how¬ 
ever, we see a systematic offset in [Ca/Fe] with respect to 


the abundances from the microlensed dwarf stars in Bensbyj 
et al. (2013]). This difference is not seen in [Mg/Fe], This off¬ 


set is most likely due to different systematics affecting stars 
in different temperature and gravity regimes. 


Acknowledgements. We are grateful for the useful comments received from an 
anonymous referee. We warmly thank the ESO Paranal Observatory staff for per¬ 
forming the observations for this programme. MZ and DM acknowledge fund¬ 
ing from the BASAL CATA through grant PFB-06, and the Chilean Ministry of 
Economy through ICM grant to the Millennium Instritute of Astrophysics. MZ 
acknowledge support by Proyecto Fondecyt Regular 1150345. DM acknowledge 
support by FONDECYT No. 1130196. 


References 

Alonso, A.. Arribas, S., & Martfnez-Roger, C. 1999, A&AS, 140, 261 
Alves-Brito, A., Melendez, J., Asplund, M„ Ramirez, I., & Yong, D. 2010, A&A, 
513, A35 

Athanassoula, E. 2005, Celestial Mechanics and Dynamical Astronomy, 91, 9 

Babusiaux, C. & Gilmore, G. 2005, MNRAS, 358, 1309 

Babusiaux, C., Gomez, A., Hill, V., et al. 2010, A&A, 519, All 

Babusiaux, C., Katz, D.. Hill, V., et al. 2014, A&A, 563, A15 

Benjamin, R. A., Churchwell, E., Babler, B. L., et al. 2005, ApJ, 630, L149 

Bensby, T., Aden, D., Melendez, J., et al. 2011, A&A, 533, A134 

Bensby, T., Feltzing, S., Johnson, J. A., et al. 2010, A&A, 512, A41 

Bensby, T„ Yee, J. C., Feltzing, S., et al. 2013, A&A, 549, A147 

Binney, J„ Gerhard, O., & Spergel, D. 1997, MNRAS, 288, 365 

Bissantz, N. & Gerhard, O. 2002, MNRAS, 330, 591 

Blitz, L. & Spergel, D. N. 1991, ApJ, 379, 631 

Bournaud, F., Elmegreen, B. G., & Elmegreen, D. M. 2007, ApJ, 670, 237 
Bournaud, F., Elmegreen, B. G., & Martig, M. 2009, ApJ, 707, LI 
Cao, L„ Mao, S„ Nataf, D„ Rattenbury, N. J., & Gould, A. 2013, MNRAS, 434, 
595 

Clarkson, W. I., Sahu, K. C., Anderson, J., et al. 2011. ApJ, 735, 37 
Combes, F., Debbasch, F., Friedli, D., & Pfenniger, D. 1990, A&A, 233, 82 
Cunha, K. & Smith, V. V. 2006, ApJ, 651, 491 
Daddi, E„ Elbaz, D„ Walter. F„ et al. 2010, ApJ, 714, LI 18 
de Vaucouleurs, G. 1964, in IAU Symposium, Vol. 20. The Galaxy and the 
Magellanic Clouds, ed. F. J. Kerr, 195 

Debattista, V. R, Carollo, C. M., Mayer, L., & Moore, B. 2005, ApJ, 628, 678 
Dekel, A. & Burkert, A. 2014, MNRAS, 438, 1870 
Dwek, E., Arendt, R. G., Hauser, M. G., et al. 1995, ApJ, 445, 716 
Elmegreen, B. G., Bournaud, F., & Elmegreen, D. M. 2008, ApJ, 688, 67 
Forster Schreiber, N. M., Genzel, R., Bouche, N., et al. 2009, ApJ, 706, 1364 
Freeman, K., Ness, M., Wylie-de-Boer, E., et al. 2013, MNRAS, 428, 3660 
Fulbright, J. P„ McWilliam, A., & Rich, R. M. 2006, ApJ, 636, 821 
Fulbright, J. P„ McWilliam, A., & Rich, R. M. 2007, ApJ, 661, 1152 
Gardner, E„ Debattista, V. P., Robin, A. C., Vasquez, S., & Zoccali, M. 2014, 
MNRAS, 438, 3275 

Genzel, R., Newman, S., Jones, T., et al. 2011, ApJ, 733, 101 
Gonzalez, O. A., Rejkuba, M., Zoccali, M., et al. 2011, A&A, 530, A54 

Gonzalez, O. A., Rejkuba, M., Zoccali, M., et al. 2013, A&A, 552, Al 10 

Gonzalez, O. A., Rejkuba, M., Zoccali, M., et al. 2012, A&A, 543, A13 

Gratton, R. G. & Sneden, C. 1990, A&A, 234, 366 
Hill, V., Lecureur, A., Gomez, A., et al. 2011. A&A, 534, A80 
Howard, C. D„ , R. M., Clarkson, W„ et al. 2009, ApJ, 702, L153 
Immeli, A., Samland, M., Gerhard, O., & Westera, P. 2004, A&A, 413, 547 
Johnson, C. I., Rich, R. M., Fulbright, J. P., Valenti, E., & McWilliam, A. 2011, 
ApJ, 732, 108 

Johnson, C. I., Rich, R. M., Kobayashi, C., Kunder, A., & Koch, A. 2014, AJ, 
148, 67 

Johnson, C. I.. Rich, R. M., Kobayashi, C., et al. 2013, ApJ, 765, 157 
Kormendy, J. & Kennicutt, Jr., R. C. 2004, ARA&A, 42, 603 
Kunder, A., Koch, A., Rich, R. M„ et al. 2012, AJ, 143, 57 
Lecureur, A., Hill, V., Zoccali, M., et al. 2007, A&A, 465, 799 
Li, Z.-Y. & Shen, J. 2012, ApJ, 757, L7 

Lindler, D. 2001, FUSE IDL Tools Reference Handbook, in Users Guide, Far 
Ultraviolet Spectroscopic Explorer 

Mancini, C., Forster Schreiber, N. M.. Renzini, A., et al. 2011, ApJ, 743, 86 

Martinez-Valpuesta, I., Shlosman, I., & Heller, C. 2006, ApJ, 637, 214 

McWilliam, A. & Rich, R. M. 1994a, ApJS, 91, 749 

McWilliam, A. & Rich, R. M. 1994b, ApJS, 91, 749 

McWilliam, A. & Zoccali, M. 2010, ApJ, 724, 1491 

Melendez, J., Asplund, M., Alves-Brito, A., et al. 2008, A&A, 484, L21 

Minniti, D. 1996, ApJ, 459, 579 

Minniti, D., Lucas, P. W., Emerson, J. P., et al. 2010, New A, 15, 433 


11 















O. A. Gonzalez et al.: The Giraffe Inner Bulge Survey (GIBS) II 


Mucciarelli, A., Pancino, E., Lovisi, L.. Ferraro, F. R., & Lapenna, E. 2013, ApJ, 
766, 78 

Nataf, D. M., Udalski, A., Gould, A.. Fouque, P., & Stanek, K. Z. 2010, ApJ, 
721, L28 

Ness, M., Freeman, K., Athanassoula, E„ et al. 2013a, MNRAS, 430, 836 
Ness, M., Freeman, K., Athanassoula, E„ et al. 2013b, MNRAS, 432, 2092 
Newman, S. F., Genzel, R„ Forster Schreiber, N. M., et al. 2013, ApJ, 767, 104 
Nomoto, K„ Thielemann, F.-K„ & Wheeler, J. C. 1984, ApJ, 279, L23 
Pasquini, L., Avila, G., Allaert, E., et al. 2000, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 4008, Optical and 
IR Telescope Instrumentation and Detectors, ed. M. Iye & A. F. Moorwood, 
129-140 

Ramirez, I. & Melendez, J. 2005, ApJ, 626, 465 

Ramirez, S. V., Stephens, A. W., Frogel, J. A„ & DePoy, D. L. 2000, AJ, 120, 
833 

Rangwala, N. & Williams, T. B. 2009, ApJ, 702, 414 

Rattenbury, N. J„ Mao, S., Sumi, T„ & Smith, M. C. 2007, MNRAS, 378, 1064 

Rich, R. M. 1988, AJ, 95, 828 

Rich, R. M. & Origlia, L. 2005, ApJ, 634, 1293 

Rich, R. M„ Origlia, L„ & Valenti, E. 2007, ApJ, 665, LI 19 

Rojas-Arriagada, A„ Recio-Blanco, A., Hill, V„ et al. 2014, A&A, 569, A103 

Ryde, N., Gustafsson, B., Edvardsson, B., et al. 2010, A&A, 509, A20 

Sadler, E. M„ Rich, R. M„ & Temdrup, D. M. 1996, AJ, 112, 171 

Saito, R. K., Zoccali, M„ McWilliam, A., et al. 2011, AJ, 142, 76 

Shen, J., Rich, R. M., Kormendy, J., et al. 2010, ApJ, 720, L72 

Smith, G. & Ruck, M. J. 2000, A&A, 356, 570 

Sneden, C. 1973, ApJ, 184, 839 

Soto, M„ Rich, R. M„ & Kuijken, K. 2007, ApJ, 665, L31 
Stanek, K. Z., Mateo, M„ Udalski, A„ et al. 1994, ApJ, 429, L73 
Stetson, P. B. & Pancino, E. 2008, PASP, 120, 1332 
Sumi, T„ Wu, X„ Udalski, A„ et al. 2004, MNRAS, 348, 1439 
Tacconi, L. J., Genzel, R„ Neri, R., et al. 2010, Nature, 463, 781 
Tacconi, L. J., Neri, R., Genzel, R., et al. 2013, ApJ, 768, 74 
Tinsley, B. M. 1979, ApJ, 229, 1046 

Udalski, A., Szymanski, M., Kubiak, M., et al. 2002, Acta Astron., 52, 217 
Valenti, E„ Zoccali, M., Renzini, A., et al. 2013, A&A, 559, A98 
Wegg, C. & Gerhard, O. 2013, MNRAS, 435, 1874 
Woosley, S. E. & Weaver, T. A. 1995, ApJS, 101. 181 
Zoccali, M., Gonzalez, O. A., Vasquez, S., et al. 2014, A&A, 562, A66 
Zoccali, M., Hill, V., Lecureur, A„ et al. 2008, A&A, 486, 177 
Zoccali, M., Lecureur, A., Barbuy, B., et al. 2006, A&A, 457, LI 
Zoccali, M., Renzini, A., Ortolani, S„ et al. 2003, A&A, 399, 931 


12 



